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The subbands of weakly coupled double-layer two dimensional electron gas systems consist of 
narrowly spaced pairs whose corresponding wavefunctions are symmetric and antisymmetric com- 
binations of isolated layer subband wavefunctions. The energetic spacing within a pair is 2t where 
t is the interlayer tunneling amplitude, t is an important parameter in modeling these systems 
and, if interactions could be neglected, it would be proportional to beating frequencies seen in 
weak-field magnetic oscillations experiments and therefore readily measurable. We point out that 
interactions alter the beating frequency. We discuss similarities and differences between this effect 
and exchange-correlation enhanced spin-splitting. 

PACS numbers: 73.20.Dx, 73.40.Hm, 75.20.En 
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The development of of high mobility double-quantum- 
well structures has opened up a fruitful area of research 
in low-dimensional physics. The many-body physics of 
these systems is especially interesting in the strong mag- 
netic field regime where quenching of the kinetic en- 
ergy due to Landau quantization leads to a variety of 
strong-correlation effectsu, including novel quantum Haj 
effectsHu and spontaneous interlayer phase coherenceai 
In models of these systems a critical parameter is the 
amplitude t for quantum tunneling of electrons between 
the layers. In a tight-binding approximation, the sub- 
bands in weakly coupled electrically balanced double- 
layer systems occur in pairs, separated in energy by 2t, 
whose wavefunctions are the symmetric and antisymmet- 
ric combinations of isolated layer subband wavefunctkins^ 
For individual samples the parameter t is routinelytluijfl 
determined by analyzing Shubnikov-de Haas oscillations 
in the resistance at weak magnetic fields. The magnetic 
oscillation frequencies are proportional to the areas en- 
closed by the circular Fermi circles for the two subbands. 
In the past it has been assumed that for double-layer 
systems the difference in cross-sectional areas of the sym- 
metric and antisymmetric subbands can be related to t 
by using the following equation which would be valid for 
non-interacting electrons: 
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Here kp and n are the Fermi wavevectors and densities 
for symmetric (S) and antisymmetric {A) subbands, m* 
is the semiconductor band structure effective mass, and 
I'Q = m */{tiTt^) is the non-interacting electron density of 
states for each two-dimensional subband. In this paper 
we point out that when electron-electron interactions are 
included Eq. |^ is no longer valid. The interlayer tunnel- 



ing parameter appearing this equation must be replaced 
by an effective value t* . The change in the effective value 
of t is analogous to the familiar exchange-correlation en- 
hancement of the spin-splitting of electronic bands in the 
presence of a magnetic field. However, as we will discuss, 
there are important differences. As a consequences of 
these, the effective splitting can be either enhanced or 
suppressed depending on the layer separation. 

We will couch our analysis in terms of a pseudospin 
representations for the layer degree of freedom which has 
proved useful in work on double-layer systems in strong 
magnetic fields. We restrict our attention here to the 
case where the two coupled layers are identical. Then 
the second quantized Hamiltonian of the system H = 
Hi + H2 , where the one body term 
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and the interaction term 
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Here s in the spin label, a is the pseudospin label, Vo{q) = 

iVA{q) + VEm/2, and VM = {VA{q)-VE{q))/2 

where Va (<z) and Ve (q) are respectively the intra and in- 
ter layer effective electron-electron interactions, (cr = ±1 
for electrons in left and fright layers and a — —cr.) 
In the narrow layer limita, which we use for illustra- 
tive calculations below, VA{q) = 2-Ke'^/q and VE{q) = 
exp (— gd) Va((?) where d is the layer separation. 

The invariance of the Hamiltonian under reversal of 
pseudospin labels guarantees that the Greens function 
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is diagonal when expressed in the representation of 
symmetric and antisymmetric states. It therefore fol- 
lows from the Luttinger theoremllSl that the relationship 
between the Fermi circle areas measured in magneto- 
oscillation experiments and the densities of symmetric 
and antisymmetric electrons is not altered by interac- 
tions. However the relationship between these densities 
and the hopping parameter which we calculate below, 
is altered. We discuss this relationship in terms of the 
linear response of the x component of the pseudospin po- 
larization to the hopping term in the Hamiltonian: 
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In this equation t" is a Pauli spin matrix. Since the hop- 
ping term can be recognized as Zeeman coupling to the 
X component of the pseudospin, the required calculation 
would be identical to that for the_apiii susceptibility of a 
two-dimensional electron systemliilO if the pseudospin- 
dependent term in H2^ which vanishes for c? — > 0, were 
not present. We will refer to as the transverse pseu- 
dospin susceptibility to distinguish it from the response 
of the z pseudospin component to a bias potential, which 
is sharply reduced by screening effects omsent already 
at the level of the Hartree approximationlj for this lon- 
gitudinal pseudospin susceptibility. We have evaluated 
the transverse pseudospin susceptibility using the gen- 
eralized random-phase-approximation (GRPA) and find 
that it is qualitatively altered by the pseudospin depen- 
dence of the interaction term in the Hamiltonian. This 
calculation is outlined in the following paragraph. The 
GRPA has important deficiencies which we discuss below, 
but it is adequate for the qualitative issues we address. 

In the GRPA the transverse pseudospin susceptibility 
is approximated by the following expression: 
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where A is the two-dimensional (2D) system area and 
Ld) is the solution to the following integral equation: 
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In Eq. ^ and Eq. |^ is the Hartree-Fock approximation 
quasiparticle energy dispersion relation at t = 0, 
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and fj: is the corresponding Fermi occupation factor. 
The starting point for the derivation of Eq. ^ is the 
fluctuation-dissipation theorem expression for the trans- 
verse pseudospin susceptibility in terms of thermody- 
namic fluctuations of the pseudospin operator: 
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The GRPA is most economically derivedEfl by making a 
Hartree-Fock factorization of the equation of motion for 
X^^ {q, t) . Alternately, the same expression can be ob- 
tained in diagrammatic perturbation theory by evaluat- 
ing ladder-sum vertex correction to the bubble skeleton 
diagram. In the limit of interest {oj — Q and q Q) the 
ratio of Fermi factors and quasiparticle energy differences 
in these equations simplifles to {m^^ /hkp)S{k — kp) 
where m^^ is the Hartree-Fock approximation quasi- 
particle effective mass allowing. With this simpliflcation 
Eq. 1^ can be solved analytically and the sum in Eq. ^ 
performed. The result is 
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Eqs. and ^ may be combined to express the 

transverse pseudospin susceptibility in a form which is 
familiar from the electron gas theory of the spin suscep- 
tibility and which form the basis of the following discus- 
sion. 
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where S is the Stoner enhancement factor and the Stoner 
interaction parameter 
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The GRPA result for / as a function of layer separa- 
tion is plotted in Fig. |^. The most remarkable feature 
of this result is that / changes sign and the susceptibil- 
ity is suppressed rather than enhanced for kpd > 1.13. 
In an electron gas, the usual spin susceptibility is in- 
variably enhanced by interactions because the band en- 
ergy cost of spin alignment is partially offset by increased 
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attractive exchange interactions between parallel spins. 
The exchange energy gain due to the alignment of trans- 
verse pseudospin components occurs when coherence is 
established between electrons in opposite layersEj and 
is weaker because it comes only from intcrlayer inter- 
actions. In the limit d oo, this vertex contribution to 
/ vanishes and we are left only with the interaction ef- 
fect on the quasiparticle dispersion in each layer. For the 
Hartree-Fock approximation vanishes in Coulom- 

bic systems, causing the GRPA value of / to approach 
—oo and S to vanish as d oo. This property is an 
artifact of the GRPA. However the fact that for d ^ oo 
S depends only on the one-particle Greens function is a 
general result. Using the property that the two layers are 
uncorrelated we have derived a formal general expression 
for this limit: 
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huj — e' + e -\- irj 

where p^{e) is the spectral weight functionQ of the iso- 
lated layer one-body Greens function normalized so that 
it integrates to unity. The quasiparticle pole contribution 
to limg--^o X^^io) for ~^ oo, which we expect to be domi- 
nant, gives a Stoner enhancement factor S = z^nfi^ /m* 
where z is the quasiparticle normalization factor, and 
rrfi^ is the quasiparticle effective mass at the Fermi en- 
ergy. For moderate density 2D electrQii.|S¥Stems this fac- 
tor is considerably smaller than unityt3E3, in qualitative 
agreement with the large d GRPA result. 

The GRPA estimate of S at general d can be improved 
by using a static screening approximation for Va{(1) a-nd 
Ve (q) ; this procedure is somewhat ad hoc but it is known 
to work reasonably well numerically for analogous quan- 
tities in single layer systems. In particular, quantum 
Monte Garlo estimatesEJ for the quasiparticle mass of 
2D electron systems are close to screened Hartree-Fock 
values. However the reduction of the ordinary spin- 
susceptibility compared to the GRPA is overestimated 
by this approximation. For q < 2kF, theJlPA screened 
intralayer and interlayer interactions arellZI 
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Results obtained for the d dependence of the Stoner fac- 
tor using these screened interactions are illustrated in 
Fig. ||. The results are only weakly dependent on density 
over this wide range which covers the density range of 
typical samples. (For 0, S* — > 1, independent of 

kpd.) Again we find that the susceptibility is suppressed 
rather than enhanced for kpd > 1.1. 



As mentioned above, many-body effects in double-layer 
quantum Hall systems are often quite sensitive to the 
precise numerical value of the interlayer hopping param- 
eter t, especially in the limit of strong magnetic fields. 
Reliable comparison between theory and experiment can 
be dependent on accurate estimates of this system pa- 
rameter. In this work we have pointed out that a non- 
interacting electron interpretation of weak field magnetic 
oscillation experiments does not necessarily yield an ac- 
curate value for t. The true value of t can be obtained 
from the apparent value by dividing by the Stoner en- 
hancement factor for the transverse pseudospin suscep- 
tibility. Over the density range of typical samples our 
results are accurately interpolated by the formula 
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where x — kpd. The precise value of S will depend on 
density, and on finite thickness effects in the 2D electron 
layers, and may differ somewhat from the GRPA static- 
screening value. 
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FIG. 1. GRPA transverse pseudospin susceptibility Stoner 
interaction parameter 7 as a function of kpd. In this ap- 
proximation likpd) / I{d = 0) is independent of density. 
^{^) ~ {'^Y^^rs/'K where Vs is the usual electron gas den- 
sity parameter defined here in terms of the density per layer, 
n = {wrsao)~^. Here ao = fi^/m*e^ is the semiconductor 
Bohr radius. For GaAs ao ~ 10 nm. The Stoner enhance- 
ment parameter S = {1 — . 

FIG. 2. GRPA transverse pseudospin susceptibility Stoner 

enhancement parameter S, calculated using statically 
screened Coulomb interactions, as a function of kpd for a se- 
ries of 2D electron gas density parameters. Results are plotted 
for r, = 0.5, 1.0, 1.5, and 2.0. 
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